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The neutron and proton single-particle spectral functions in asymmetric nuclear 



^ ' matter fulfill energy weighted sum rules. The validity of these sum rules within the 

lO ' self-consistent Green's function approach is investigated. The various contributions 

o _ 

\ to these sum rules and their convergence as a function of energy provide information 



about correlations induced by the realistic interaction between the nucleons. These 
features are studied as a function of the asymmetry of nuclear matter. 

PACS numbers: 21.65.+f, 21.30.Fe 

I. INTRODUCTION 

The microscopic description of the single-particle properties in nuclear matter must deal 



with the treatment of nucleon-nucleon (NN) correlations [Ij, |2|. In fact, the strong short- 
range and tensor components which are required in realistic NN interactions to fit the NN 
scattering data, lead to corresponding correlations in the nuclear wave function. These NN 
correlations are important to describe bulk properties of dense matter. They also modify the 
spectral distribution of the single-particle strength in a significant way. Recent calculations 
have demonstrated without ambiguity that NN correlations produce a partial occupation of 
the single-particle states, which are completely occupied in the mean field approach, and a 
wide distribution of the single-particle strength in energy. These two features have found 
experimental grounds in the analysis of (e, e'p) reactions j^Q]. 
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Historically, various tools have been employed to account for correlations in the nuclear 
many-body wave function. These include the traditional Brueckner hole-line expansion 
and variational approaches using correlated basis functions 0]. Attempts have also been 
made to employ the technique of a self-consistent evaluation of Green's functions (SCGF) 
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to the solution of the nuclear many-body problem. This method offers various 
advantages: (i) the single-particle Greens function contains detailed information about the 
spectral function, i.e., the distribution of single-particle strength as a function of missing 
energy and momentum; (ii) the method can be extended to finite temperatures, a feature 
which is of interest for the study of the nuclear properties in astrophysical environments; 
(iii) the Brueckner Hartree-Fock (BHF) approximation, the approximation to the hole-line 
expansion which is commonly used, can be considered as a specific prescription within this 
scheme. 

Enormous progress in the SCGF applications to nuclear matter has been achieved in the 
>a.t yea., botK at S and fin.te te.p.atu. H y Q Q. The efforts at T ^ 

are mainly oriented to provide the appropriate theoretical support for the interpretation of 
the (e, e'p) experiments, while research at finite T is essentially focused on the description 
of the nuclear medium in astrophysical environments or on the interpretation of heavy ion 
collisions dynamics. In all key quantity is the single-particle spectral function which 

measures the possibility to add or remove a particle with a given momentum at a specific 
energy. A useful way to study the properties of these single-particle spectral functions is 
by means of the energy weighted sum rules, which are well established in the literature and 
which have been numerically analyzed for symmetric nuclear matter at zero ^ and finite 
temperature jisl . 

Recently, the single-particle spectral functions in hot asymmetric nuclear matter have been 
calculated within the SCGF framework In this computation, the bulk properties of 

asymmetric dense matter were computed, namely the energy per particle, the symmetry 
energy and the chemical potentials. In addition, and since the SCGF approach gives access 
to the correlated momentum distributions of neutrons and protons, the influence of the 
asymmetry on the depletion of the momentum states could be studied. In neutron-rich 
matter the proton states with momenta below the Fermi momentum are more strongly 
depleted than in the symmetric system at the same density. In this case the occupation of 
neutron states within their Fermi sphere are enhanced compared to the symmetric case. This 
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indicates that the proton-neutron interaction is a stronger source of correlation as compared 
to proton-proton and neutron-neutron interaction. 

In this paper, we want to use the energy weighted sum rules to investigate the dependence 
of the single-particle spectral functions on the asymmetry more in detail. These sum rules 
may help to explore the isospin dependence of the short-range correlations in asymmetric 
matter. Besides, the analysis of the energy weighted sum rules provides valuable tests on 
the numerical accuracy of our calculations. All the computations discussed in this paper are 
performed in the framework of SCGF employing a fully self-consistent ladder approximation 
in which the complete spectral function has been used to describe the intermediate states 
in the Galistkii-Feynman equation 12. i^. 

After this introduction, we outline the derivation of the energy weighted sum rules in Sec- 
tion 2. The results obtained for hot asymmetric nuclear using the charge dependent Bonn 
potential CDBONN Q are presented in Section 3, where a short summary of the main 
conclusions is given as well. 



II. SUM RULES FOR ASYMMETRIC NUCLEAR MATTER 

For a given Hamiltonian if, the single-particle Green's function for a system at finite tem- 
perature can be defined in a grand-canonical formulation: 



ig. 



-(kt; k't') = Tr {pT [a^^{t)ai,^{t')\ } , (1) 



where the subindex t = n{p) refers either to neutrons or protons and T is the time ordering 
operator that acts on a product of Heisenberg field operators akr(^) = e**^akre~**^, in such 
a way that the field operator with the largest time argument is put to the left. The trace is 
to be taken over all energy eigenstates and all particle number eigenstates of the many-body 
system, weighted by the statistical operator, 

p = ]_^-P{H~tj.„N„~fipNp) _ ^2) 

Z 

(5 denotes the inverse of the temperature and /i„ and /ip are the neutron and proton chemical 
potentials associated to a given average number of neutrons and protons. A^^^ and Np are 
the operators that count the total number of neutrons and protons in the system, 

Nr = Y.<^ir{t)a^r{t). (3) 
k 



4 



The normalization factor in Eq. Q is the grand partition function of statistical mechanics, 

For a homogeneous system, the Green's function is diagonal in momentum space and depends 
only on the absolute value of k and on the difference ^ = t' — t. Starting from the definition 
of the Green's function, we first focus on the case ^ > 0. In order to recover the expression 
in the grand canonical ensemble average of the occupation number nr{k) for C, = 0~^, the 
following definition of the correlation functions includes an additional factor —i with 
respect to the definitions of the Green's function g^, in Eq. (jT)), 

g<{k, = Tr {pe^«^aLe-^«^ak.} . (5) 

Qui^^O be expressed as a Fourier integral over the frequencies, 

_ (6) 

where A'^{k^uo) is defined by 

A<{k,uj)=2T:Y. I I auu I \^ 5{uj - {Er - Es)) . (7) 

A similar analysis can be conducted for ^ < 0, yielding a function 

^>(fc,c^) = e'^("-'^'')A<(A;,cu). (8) 

Finally, the single-particle spectral function is defined as the sum of the two positive func- 
tions, A"^ and A^ , 

A,{k,u:) = A<{k,u) + A>{k,uj) . (9) 

In our computations, we obtain the full spectral function Ay{k,uj) and we make use of the 
relations: 

A<{k,uj) = f{uj)A,{k,u;) , A>{k,uj) = (1 - f{u))A,{k,uj) (10) 

to obtain and A-^ , where f{uj) is the Fermi-Dirac distribution function. The integration 
of the function A"^ leads to the occupation probability 

n.{k)= r^^A<ik,co), (11) 

for the state with isospin u and momentum k at the inverse temperature /?. 
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The functions and can be compared with the hole and particle part of the spectral 
function at zero temperature. In particular, the limit of A"^ at zero temperature reads 

4 = 2vrE I (*^^ I a,. I ^0 ) I' ^(^ - (^^o^ - Ef-')) , (12) 

r 

where | "^q) is the ground state with Nn neutrons and Np protons (such that A = Np + Nn) 
and I '^^~^) labels the excited energy state of a system with a neutron or a proton less 
(depending on which type of annihilation operator aku has been applied to the ground 
state). By its own definition, it is clear that the lowest possible energy of the final state is 
the ground state energy of the A — 1 particle system, so that there is an upper limit for the 
hole spectral function uj^ = Eq — Eq"^" = /i^, both for neutrons and protons. In a similar 
way one can define the particle part of the spectral function and find a lower bound for the 
excitation energy of the A + 1 particle system, measured with respect to the ground state of 
the A particle system. At zero temperature, the existence of these lower and upper bounds 
causes a complete separation in energy between the particle and the hole part of the spectral 
function. 

In the T matrix approximation to the self-energy reported in Ref . , one can determine the 
single-particle Green's function as the solution of Dyson's equation for any complex value 
of the frequency variable z 

The analytical properties of the finite temperature Green's function allows one to derive the 
corresponding Lehmann representation, which for slightly complex values of the frequency 
can be written as 

g,{k,uj\iri)= / • (14) 

J-oo lT\ IjJ — 17] 

The sum rules for the spectral functions can be derived from the asymptotic behavior at 
large uj by expanding the real part of both previous expressions for the Green's function, 
Eqs. (fTSj) and (HH), in powers of Xjuj. This yields 



Re gv{k, uj) = —\ \ ^ — 

U UJ 



h lim ReT,^(k, oj 

2m 



+ ... (15) 



and 



Re g^{k, u) = - \ dJA^{k, uj') + - du'u'A^{k, + ... \ . (16) 
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By comparing the first two expansion coefficients, one finds the mo 

/■+°° duj 

/ ^A^{k,uj) = l, (17) 

and the mi sum rules 

— u;A^(A;,cj) = h hm ReSj,(A;, c^;) . (18) 

Is is worth to mention that in our scheme the self-energy is derived in the T matrix ap- 
proximation and so its real part is computed from the imaginary part using the following 
dispersion relation: 

ReE.(^.^)^SrW--r°°dA '°'^-''-V"'' . (19) 

71 J~oo UJ — A 

In the derivation of the previous equation, the spectral decomposition of the Green's function 
is already used, so it is a property of the T matrix approach that it fullfils the sum rule. 
Nevertheless, the sum rules still provide a useful consistency check for the numerics. The 
first term on the right hand side of Eq. (fT^ is the energy independent part of the self-energy, 

^^(^) = E / 1^ (k^kVI V I Wr)^ n.(k') (20) 

which can be identified with the limit lim^^oo ReSj,(/c, cj), since the dispersive part decays 
like \/uo for uj ±oo . Eq. ()20|) looks like a Hartree-Fock (HF) potential. However, 
Uyik) is the momentum distribution of Eq. (fTTj) containing the depletion effects due to NN 
correlations and temperature. 



III. RESULTS AND DISCUSSION 

All the results discussed in this paper have been obtained with the charge-dependent Bonn 
(CDBONN) potential. Since we want to focus our study on the dependence of the observables 
on asymmetry, we will consider only one density p = 0.16 fm~^ and one single temperature 
T = 5 MeV. This temperature is low enough to allow for conclusions in the limit of zero 
temperature but high enough to avoid instabilities associated with neutron-proton pairing 
[itI lisl. Iioj l. As for the proton fractions, we will consider three different cases. The first 
one corresponds to symmetric nuclear matter Xp = Pp/p = 0.50 and serves as a guide-line 
for the other cases. The last one has a very low proton fraction Xp = 0.04 and corresponds 
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to the /9— stable composition of matter with nucleons, electrons and muons at T = 5 MeV 
and p = 0.16 fm~^. This composition can be computed thanks to the fact that we know 
the asymmetry dependence of both the neutron and the proton chemical potentials at this 
density Finally, we consider an intermediate fraction Xp = 0.30 which is useful to 

identify the effects for low asymmetries of the system. 

We start by discussing the momentum dependence of the single-particle spectral functions 
of neutrons and protons. Depending on how far above or below the considered momentum is 
from the Fermi momentum of the corresponding particle, we expect a very different behavior 
for the spectral function. Fig. ^ shows the neutron (left panels) and proton (right panels) 
single-particle spectral functions for three different momenta {k = 0, kp and 2kp, with kp 
the Fermi momentum of each nucleon species) at a proton fraction of Xp = 0.04. The dotted 
vertical line corresponds to u = uj — fi,^ = —^y and indicates the point below which the 
contribution to the sum rule mi in Eq. (|18p is negative. 

The contributions and are given by the dashed and dash-dotted lines, respectively. 
Since we are considering finite temperature, thermally excited states are always included in 
the definition of the spectral function Eq. (0). This leads to contributions to at energies 
uj larger than /i,^ {uj larger than zero). Similarly, A^{k,uj) extends to the region below p^. 
In general there is no longer a clear separation in energy between A^ and A^, as it is the 
case for T = 0. Actually, the maxima of A^ and A^ can even coincide. 
For the case of neutrons at A; = 0, the peak of the spectral function is provided by A^. This 
is due to the fact that the position of this peak, which can be identified with the quasi- 
particle energy e^pik), is well below the chemical potential which implies that the value 
of Fermi-Dirac distribution function /(cu) at this energy oj = e^p(O) is very close to one (see 
Eq. irUll. Since the quasi-particle peak is well below the thermal effects do not fill up the 
minimum in the spectral function at = and we still observe some kind of separation in 
energy between and A^. However, around u; = there is an energy interval, in which 
both A^ and A^ are small but different from zero, i.e. an interval in which these functions 
overlap. 

A similar situation is observed for the neutron spectral function at A; = 2kp. In this case 
the quasi-particle energy is well above fin, which means that (1 — f{ui)) is very close to one 
at this energy and the peak structure is supplied by A^ (see Eq. ITIH) . Therefore we observe 
a clear separation between the hole part (A^) and the particle part (A^) of the spectral 
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function even at finite temperature. 

At k = kp, the quasi-particle peak is close to the chemical potential. This means that 
f{€qp{k)) is around 0.5 and therefore and have a peak at u; = 0. The thermal effects 
fill up the zero that the spectral function has at = in the T = limit and the overlap 
region of A^ and A^ is enhanced. 

In the case of very small proton fraction, like Xp = 0.04, the Fermi momentum kp is rather 
small and consequently the quasi-particle energies are close to the Fermi energy fip for 
k = k^ a.s well as = and k = 2kp. Therefore, although the relative distances to the 
Fermi surfaces are the same, there are strong overlaps between Ap and Ap at all the three 
momenta. Even in the case oi k = a small peak structure for Ap around u) = can be 
observed. 

In Table m (neutrons) and Table |n] (protons) we report the fraction of the integrated strength 
of below and above the corresponding chemical potential, /ijy. In fact, if we identify A^ 
with the T = hole spectral function, the integrated strength above the chemical potential 
would be exactly zero (since A"^ = for > fi^). In this sense, the integrated strength 
above the chemical potential can be considered as a genuine thermal effect. As expected, 
these thermal effects are important around kp, where the overlap between A"^ and A^ is 
significant. Notice also that for this very neutron-rich system, neutrons are less affected by 
temperature, while protons (which have a substantially lower Fermi momentum and can thus 
be considered as a dilute system) are much more influenced by temperature. For instance, 
protons show a large amount of strength above /i^ up to momenta k/k^ ~ 2. For the sake 
of completeness, the corresponding occupation numbers for each species are also listed in 
those tables. 

The mo sum rule is fulfilled with an accuracy better than .1 % for both neutrons and 
protons in the complete momentum range. Results for mi both for neutrons (upper panels) 
and protons (lower panels) are reported in Fig. |2lat the three proton fractions Xp = 0.5,0.3 
and 0.04. The vertical dotted lines indicate the location of kp for each case. Since the sum 
rule is satisfied better than 1 % in all cases, the left- and the right-hand sides of Eq. ()18|1 
(solid lines) lie on top of each other and cannot be distinguished. To understand the k- 
dependence of the different contributions to mi, it is useful to keep in mind the proton and 
neutron chemical potentials at the different concentrations given in Table ITTTl 
The lower dash-dotted line shows the mi contribution from A"^. For momenta below kp, the 
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contribution of A"^ to mi is dominated by the quasi-particle peak, which hes below fi^. As 
the momentum increases, the peak appears closer to = and its weight in the integral is 
diminished, which makes the integral smaller in magnitude and thus the contribution to the 
sum rule becomes an increasing function of k. When we get closer to the Fermi momentum 
kp, the peak moves to u fi^ and the position of the chemical potential becomes crucial. 
As far as the chemical potential is negative, the contribution of A"^ will be negative because 
the Fermi-Dirac function falls off close to fi^, and the integral will only catch a little positive 
zone if the chemical potential is close to ~ 0. If, on the other hand, the chemical potential 
is positive, the integrand is not zero for > and there can be a non-negligible positive 
contribution to the integral when the quasi-particle peak lies between u = and uj = fi^. 
This cancellations between positive and negative contributions give rise to the structures 
observed for the sum rule close to the Fermi momentum. In the upper-right panel of 
Fig. 121 for instance, the asymmetry is so extreme that the neutron chemical potential > 
and the positive contribution to the integral is important enough to pull the sum rule to 
positive values for k ~ kp. Finally, for high momenta, A"^ is strongly suppressed (the quasi- 
particle peak lies in the region u > fi^, which is suppressed by the Fermi-Dirac factor) and 
the contribution to the sum rule goes to zero. 

The upper dash-dotted line displays the contribution from A^ . Due to the short-range cor- 
relations, there is always a high energy tail that gives rise to a positive contribution which, 
for momenta well below kp, is nearly constant. In the case of the neutron sum rules, this 
constant decreases slightly with decreasing proton fraction, although the chemical potential 
fin increases significantly. Such a small decrease is in accordance with the fact that the 
neutron-proton correlations are dominant compared to the neutron-neutron ones. These 
high energy tails (which are mainly caused by short-range correlations) are, nevertheless, 
not much affected by asymmetry. When the momenta becomes larger than kp, mi increases 
monotonously following the location of the quasi-particle peak, which moves to higher en- 
ergies when the momentum grows. 

Also for protons one observes a contribution from A-^ to the sum rule which is almost a 
constant for momenta smaller than the Fermi momentum. In this case, however, this con- 
stant is almost independent of the proton fraction, although the chemical potential fip gets 
significantly more attractive with decreasing Xp. At small values of Xp, the neutron abun- 
dance is large and this leads to strong correlation effects in the proton spectral function. In 
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particular at Xp = 0.04, which is the case corresponding to the lowest Fermi momenta, the 
total contribution of is the result of a balance between positive and negative contribu- 
tions. These can result in a minimum for k ~ kp. In that particular proton fraction (see 
Fig. P), the quasi-particle peaks lies always below u = and thus it gives rise to a negative 
contribution. The relative width of this peak and the detailed structure of the high-energy 
tail can push the integral to negative values (almost —15 MeV at ~ 160 MeV). At larger 
momenta, however, the contribution to mi of starts to grow steadily. Notice that for 
this asymmetry we can not say that this is due to the movement of the quasi-particle peak, 
because even at = 2k^ the quasi-particle peak lies in the region where the contribution 
to the integral is negative. The growth in k is, to a large extent, a consequence of the 
contributions to the spectral functions at large energies. 

In addition to these contributions, we have plotted in the dotted lines of Fig. |21the Hartree- 
Fock approximation of mi at the same temperature, density and proton fraction. This 
approximation turns out to give a very good estimate of mi, an interesting result which has 



been already observed in the previous analysis for symmetric nuclear matter [IJ, |1^ . This 
result allows for a quantitative estimate of the amount of correlations produced by a given 
NN potential without performing sophisticated many-body calculations. 
Furthermore, it is worth noticing that the value of mi (specially for large k) is not so much 
affected by the asymmetry of the system. It is very similar for both neutrons and protons 
in all the three cases considered. This is a somehow surprising result, since we have seen 
that the separate contributions of A"^ and A^ can indeed be very different. Of course, at 
high momentum the kinetic energy is dominant and this could explain this fact in part. 
However, even at zero momentum the effect of the asymmetry is rather moderate, as we 
can see from Fig. El where S^, (the Hartree-Fock approximation to the self-energy) 
and the quasi-particle energy e'^^ at zero momentum are plotted as a function of the proton 
fraction. Let us consider the highly asymmetric case of Xp = 0.04. The isospin splitting 
of the Hartree-Fock energies [k = 0) at this asymmetry is around 10 MeV, which is rather 
small compared to the isospin splitting of the quasi-particle energies, which is around 55 
MeV. This means that the stronger binding of the proton states as compared to the neutron 
states is to a small extent due to the attractive neutron-proton interaction in the bare NN 
potential, the Born term of the T-matrix. The obtained attraction is then mainly due to the 
terms in the T-matrix coming from the second and higher orders in the NN potential. In 
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other words, it is an effect of strong correlations in the isospin equal to zero NN channels (like 
the ^5*1 —^Di partial wave), which lead to the deeply bound quasi-particle energy for protons 
at this large neutron abundance. These correlation effects, however, are also important for 
a redistribution of single-particle strength for protons with /c = to energies well above the 
Fermi energy, leading to a low occupation probability (see Table and a large positive 
contribution to the sum rule mi. This positive contribution yields to a value of which 
is even above the HF energy. The effects of correlations for the corresponding energies of 
the neutron state are much weaker, indicating that correlation effects in the isospin one NN 
channels are weaker. 

So we observe a compensation of correlation effects in the energy weighted sum rule. On the 
one hand, correlations lead to a more attractive quasi-particle energy (compared to the HF 
result) and therefore to a more attractive contribution of the quasi-particle pole to the sum 
rule. On the other hand, correlations shift single-particle strength to high energies, which 
yield a repulsive contribution to mi. This compensation of correlation effects explains that 
the value of the energy integrated sum rule for the correlated system can be very similar to 
the HF result. 

This compensation of correlation effects can also be observed in Fig. HI which shows the 
exhaustion of the sum rules tuq (left panel) and mi (right panel) at a momentum k = 400 
MeV for both neutrons (solid lines) and protons (dashed lines), at the three different proton 
fractions. The first thing to be observed is that the asymmetry affects mainly the amount of 
strength exhausted at intermediate energy regions, while there is not a noticeable difference 
at high energies. For Xp = 0.5, i.e. the symmetric case, the quasi-particle energies for 
neutrons and protons aX k = 400 MeV, are the same: e^p = 36.45 MeV. When the proton 
fraction decreases, the quasi-particle energies split. Neutrons become more repulsive (e^^ = 
41.88 MeV), in contrast with protons (e^^ = 30.86 MeV). This behavior is confirmed at 
Xp = 0.04, resulting in e^^ = 49.35 MeV and e^p = 25.22 MeV. In the case of mo, since the 
quasi-particle peak of protons is located at lower energies, the amount of strength exhausted 
at lower energies is larger for protons than for neutrons. This trend changes at intermediate 
energies and both strengths merge together for large energies, which is an indication of the 
fact that, for these momenta and asymmetries, the quasi-particle contributions to the sum 
rule of neutrons are larger than those of protons, thus explaining why the increase in mo at 
the quasi-particle pole is larger in the former case. The same type of analysis is valid for 
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mi. Notice that, as previously mentioned, the final value of mi is not much affected by the 
asymmetry and the behavior at high energies is the same for all the cases. 
To summarize, we have analyzed the behavior of the energy weighted sum rules of single- 
particle spectral functions of hot asymmetric nuclear matter. The sum rules are very well 
fulfilled, because the T-matrix approximation itself respects the analytical properties of both 
the self-energy and the Green's function, in which the sum rules are based. Nevertheless, 
they are a good test of the numerical consistency of the calculation, which may be helpful 
e.g. in deciding the best distribution of the energy mesh points when one has to work with 
spectral functions. 

Employing realistic NN interactions, an important source of correlations come from the 
strong components in the neutron-proton interaction. As a consequence, one observes in 
neutron-rich matter a larger depletion of the occupation probabilities for protons with mo- 
menta below the Fermi-momentum than for neutrons. One also finds that these correlations, 
at this asymmetry, lead to much more attractive quasi-particle energies for protons than 
those obtained in the HF approximation. The same effect is observed in neutrons, although 
it is considerably weaker. This shift of the quasi-particle energies tends to lead to more 
attractive energy weighted sum rules mi. This effect is compensated by the fact that the 
very same correlations are also responsible for a shift of single-particle strength to high pos- 
itive energies. As a consequence, the energy integrated sum rules and the isospin splitting 
in asymmetric nuclear matter for the correlated system yield results which are very close to 
the HF ones. In both cases, however, there is not a strong dependence on asymmetry due 
to the lack of tensor correlations. This relocation of single-particle strength, on the other 
hand, can be nicely observed in the convergence of the energy weighted sum rule, which is 
therefore an indicator of correlation effects. 
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1 / 7 17, 


beiow [ 


7oJ above /i^ 


7oJ n{k) 


0.0 


99.98 


0.02 


0.971 


0.5 


99.96 


0.04 


0.960 


1.0 


55.58 


44.42 


0.464 


1.5 


98.48 


1.52 


0.006 


2.0 


99.76 


0.24 


0.001 



TABLE I: Strength distribution of ^< aX. p = 0.16 fm"^, T = 5 MeV and proton fraction Xp = 0.04. 
The numbers give the fraction of the integrated strength above and below the neutron chemical 
potential Hn- The last column reports the occupation of the respective neutron momentum state. 



h / kFp 


below pp [ 


%] above /ip [ 


%] n{k) 


0.0 


97.62 


2.38 


0.605 


0.5 


94.88 


5.12 


0.531 


1.0 


21.62 


78.38 


0.280 


1.5 


24.26 


75.74 


0.056 


2.0 


64.10 


35.90 


0.009 



TABLE II: Strength distribution of at p = 0.16 fm~'^ T = 5 MeV and proton fraction Xp = 0.04. 
The numbers give the fraction of the integrated strength above and below the proton chemical 
potential /ip. The last column reports the occupation of the respective proton momentum state. 



Xp 


^ln [MeV] /ip [MeV] 


0.50 


-23.46 


-23.46 


0.30 


-2.78 


-48.33 


0.04 


19.21 


-91.45 



TABLE III: Neutron and proton chemical potentials at p = 0.16 fm ^, T = 5 MeV and different 
proton fractions. 
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FIG. 1: (Color online) Neutron (left panels) and proton (right panels) single-particle spectral 
functions at p = 0.16 fm~^, T = 5 MeV and proton fraction Xp = 0.04 for three different momenta. 
Both (dashed lines) and (dot-dashed lines) are displayed together with the total spectral 
function Ay (solid lines). 




FIG. 2: (Color online) Energy weighted sum rule mi (solid lines) for neutrons (panels above) and 
protons (panels below) at p = 0.16 fm""', T = 5 MeV and several proton fractions. Both the right- 
and the left-hand sides of Eq. ()18|) are displayed, but the sum rule is so well fulfilled that they 
cannot be distinguished. The contributions to nii that come from and are indicated by the 
upper and the lower dash-dotted lines, respectively. The result for the energy weighted sum rule 
obtained from the Hartree-Fock approximation is represented by the dashed lines. 
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FIG. 4: (Color online) Saturation of the sum rules mo (right panel) and mi (left panel) for neutrons 
(solid lines) and protons (dashed lines) at three different proton fractions. The momentum is k=400 
MeV and the density and temperature, the same same as in the previous figures. 



